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Abstract 

We study to one-loop order the renormalization of QCD in the Coulomb gauge 
using the Hamitonian formalism. Divergences occur which might require counter- 
terms outside the Hamiltonian formalism, but they can be cancelled by a redef- 
inition of the Yang-Mills electric field. 
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1 Introduction 

We study the renormalization of QCD in the Coulomb gauge Hamiltonian for- 
malism. By Hamiltonian form, we mean that the Lagrangian contains only 
first order terms in time derivatives, and depends upon the conjugate momen- 
tum field Ef as well as the (transverse) gluon field Af (here a is the colour 
index and i = 1, 2, 3 is a 3- vector index). This form has a number of attractive 
features: 

(i) As a Hamiltonian exists, the theory is explicitly unitary, without the 
necessity to cancel unphysical degrees of freedom with ghosts. 

(ii) The Lagrangian form of the Coulomb gauge has "energy divergences" in 
some of its Feynman integrals, that is integrals of the form (we use K for the 
spatial part of the 4-vector k) 



where / does not decrease as ko — ► oo (for fixed K). These divergences cancel 
between different Feynman graphs |Tj , but this cancellation has to be organized 
"by hand" . In the Hamiltonian form, each individual Feynman graph is free of 
such divergence. Formally 'energy divergent' integrals such as 



are assigned the value zero. 

(iii) It has been argued [2] that the Coulomb gauge throws light on con- 
finement. Certainly it is known [3] that, in the Coulomb gauge, the source of 
asymptotic freedom lies in the Coulomb potential. 

In spite of (i) above, to 2-loop order, mild energy-divergences remain [4], [5], 
[6] which result in ambiguities which have to be resolved by a prescription. This 
is connected with questions of operator ordering [7] . 

For other applications of the Coulomb gauge, for example to lattice QCD, 
see [S], 0. 

The question addressed here is the following. Ultra-violet divergences exist 
which seem to require the existence of counter-terms containing second order 
terms in time derivatives, (dAf /dt) 2 . Do these take us out of the Hamiltonian 
form? We argue that this does not happen because the divergences concerned 
can be cancelled by a redefinition of the field. 

We do not use quite the strict Hamiltonian formalism. We retain the auxil- 
liary field Aq, which contains no time derivatives and should be integrated out 
to give a nonlocal Coulomb potential term in the real Hamiltonian. It seems 
to be convenient, for the purposes of renormalization, to retain Aq in the La- 
grangian. Because of this, there is a ghost field, but it has an instantaneous 
propagator, and so is not relevant to unitarity. Its purpose is only to cancel out 
closed loops in the Aq field. 




(1) 
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2 The Feynman rules 

The Lagrangian for the Coulomb gauge is 



C' = C-±-(d % AlY (3) 



£ — • Fij — -(Ei) 2 + Ej • F 0j 



(where a will eventually tend to zero to go to the Coulomb gauge), 

—F^ ■ F^ — —{ 

+diC*diC + gdiC* ■ (A, A c) 
+Ui ■ [die + g(Ai A c)] 
+u • [d c + g(A A c)] 

-^K-(cAc)+ S v r (E,Ac) (4) 

where we use a colour vector notation, and 

F?j=d i A a j-djA a i +gf abc A b i A C j 

and 

(A; A c) a = f abd A\c d (5) 

Here c,c* are the ghost fields, and the sources Uj,v n and K are inserted for 
future use in formulating the BRST identities. The conjugate momentum (elec- 
tric) field E m could be integrated out to obtain the ordinary Lagrangian for- 
malism, but for the Hamiltonian formalism it must be retained. 

We will use indices m,n, ... — 1,2,3 to denote the (spatial) components of 
E, so the seven fields are (Af, Aq, E%). We will use indices /, J, .. to denote the 
seven indices (i,0, n). The bilinear part of the Lagrangian in momentum space 
is a 7 x 7 matrix 

( -K\Tij+Lij/a) 
Su5 ab = 

V ikoSm-j -iK m 

where 

Tij ~ Sij Lij , Lij = KiKj I K , 

fc 2 = fc 2 - K 2 . (6) 
For the propagators, we need the inverse 

/ T t j / k 2 — aL t j / K 2 ak K t /{K 2 ) 2 -ik T m /k 2 \ 

S7j 1 Sab=\ akoKj/iK 2 ) 2 l/K 2 + ak 2 /(K 2 ) 2 iK n /K 2 . (7) 

V ik T mj /k 2 -iK m /K 2 T mn K 2 /k 2 J 

We can now let a — * 0, to obtain the Coulomb gauge. From this, and the 
interaction terms in the Lagrangian (4), we can read off the Feynman rules. We 
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represent the Aj field by dashed lines, the E„ field by continuous lines, and the 
A field by dotted lines. With this notation, we now list the rules (a factor 
of jrfppq is to be included for each propagator, and a factor of (27r) 4 i for each 
vertex). If we choose the propagators in fig.l to be the negative of the matrix 
(7), the extra factors of ^) 4 i ^ or ^ nc propagator and (27r) 4 i for the vertices 
cancel. 



3 The ultra-violet divergences 

The divergent graphs with 2 and with 3 external lines are shown in Figures 4 
till 31. Examples of the method of evaluation of divergent parts are given in 
Appendices A and B. 

The ultra-violet divergent parts of these graphs are, in terms of the divergent 
constant (using dimensional regularization in 4 — e dimensions) 

C G T(e/2), (8) 
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(where the superfix (4), (5) etc. refers to the corresponding figure and ITjj, 
n i...n mn denote self-energies, Vijk, Vom-.-Vom vertices and A stands for dia- 
grams with external ghost lines), are: 

Ilg )Qb = z C [Ifc % + - KiKj\5 ab (9) 



L 3 

n-Q )afc = -^ick KiS ab 

n(6)ab 1 . (/ 2j 
00 = 3 tcK °ab 

n (7)afc = 
nil " 

n (9)afc _lv c A X K 
LS -mn g mn a " 



vS k ° )abc (p, q ,r) = 


-\cgf abc [{Q 


- P) k S tJ + (R 


- Q)iS jk + (P 


- R)jS lk 


vS k 1)abc ( P , q ,r) = 


-^gf abc [{Q 


- P) k 6 ij + (R 


- Q)iS jk + (P 


— R)jSik 


Vg 2)abc (p, q ,r) = 


-\cgf abc [(Q 


- P) k 8 ij + (R 


- Q)iS jk + (P 


- R)j5ik 


V^ )abC (p, q ,r) = 


\cgf abc [{Q- 


- P) k Sij + (R - 


- Q)iS jk + (P - 


- R)j$ik] 



V^ )abc { P , q ,r) = \cgf^{R-Q) t 



10) 

11) 
12) 
13) 

14) 

15) 

16) 

17) 

18) 

19) 
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t& 5)o6c (p, g, r) - -|c 5 / o6c (i2 - Q), (20) 



Voo (P> r ) = o c 9f abc (R ~ Q)% (21) 



'too IPiQ'iH — g 

r (17)a6c, „ „n _ 1 „„fabc( 



V^' aoc (p, q, r) = -cgf abc (R - Q)i (22) 

^^(p.g.r^O (23) 

^ 9)a6c (P,9,r-) = ^c ff r bc (i?-Q) (24) 

<f o6c (p, 9, r) = -\cgf abc {R - Q) (25) 

V$ 1)o6c (p,«,r)=0. (26) 
Graphs involving external E TO line are 

V™ abc (p,q,r) = ±icgf° b °6 im (27) 

^°o )a6c (p ) g,r) = -^c 9 rt ro (28) 

V^ o6c (p,«,r)=0. (29) 

All other graphs involving external E m -lines are convergent. The divergent 
parts of graphs with open ghost line are 

A {22)ab (q) = -\icQ 2 5 ab (30) 

Af 3)Qb (g) = -^cQ j( 5 ah (31) 

A^ abc (p,q) = (32) 

A^ abc (p,q,r) = (33) 

A 26)abC M=C. (34) 

Af )abc (p,q)=0 (35) 

Al 28 ) abc (p,g) = 0. (36) 
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4 Counter-terms 

Let 



xL{x) 



be the original action, T be the complete effective action, and let T\ be 
effective action to one-loop order. The complete BRST identities are 

^ „ dT dT or or or dT or or n 
r * r = 1 1 — 1 — = o. 

dAi dui dA <9u 9c dK d~Ei dwi 
So to one-loop order 

Ti *r + r *Ti = Ari = o 



where 



dT d dT d dT d dT d 



dAi dui dui dAi dA a du du <9A 
dT d dT d dT d dT d 



and 



A 2 =0. 

One class of solutions to this equation is of the form 

if > = AG, 

where the allowed form of G is, in terms of constants a 5 , ...an, 

G = a 5 Ai ■ (uj + d l c*) + a e A ■ u + a 7 c • K + a 8 Ei • Vj 

+a 9 Vi • di A + aioVi • d Ai + anVj • (A A A;). 
Other solutions of equation (39) are the explicitly gauge-invariant terms 

if = ai(F 4 ,) 2 + a 2 E 4 • F M + a 3 (F 0t ) 2 + a 4 (E,) 2 . 

Finally, by differentiating (38) with respect to the coupling constant g 
specialising to one-loop order, we see that 

Ar («i) = 

where (do being another divergent constant) 

„(»*») dT 

Combining these three contributions, we obtain 

ri - if > + if iJ + rf 4) - J d A x c(x) 



where 



£ 1 = ai(F ij ) + (a 2 + a 8 + a 9 )E l ■ F M 
+(a 3 - a 9 )(F 0l ) 2 + (a 4 - a 8 )(E l ) 2 

+a 5 Fij ■ djAi - (a 5 + -a )gF i:j ■ (A, A Aj) 

-(a + a 5 + a 6 )5E l • (Aj A A ) + Ej • (a 5 <9 Aj - a 6 d t A ) 
-a 5 (uj + <9jC*) • d l c + a gd l c* ■ (A, A c) 
-a 6 u • 9 c + a gu • (A A c) 
-a 7 (ui + diC*) ■ {diC + g(Ai A c)} 
+a gui ■ (Aj Ac) - a 7 u • {<9 c + .g(A Ac)} 

+ \g{a-7 - «o)K • (c A c) + (a - a 7 )gvi ■ (Ej A c). (48) 

The conditions coming from the vanishing ghost graphs Figs. 24, 25, 26, 27 
and 28 are particularly simple. They fix 

ctg = — ctio 

an = -ga$ 

a Q = a 7 = -a 6 . (49) 

In order for the counter-terms to cancel the divergences in the other graphs, we 
require the conditions 

4ai — 2<X5 = — c 
4ai - 3a 5 - a = 

1 

6 
4 

«6 - 0,5 = gC 

4 

a-5 + a 7 = --c 
2 

a 4 - a 8 = -c 

a-2 + a 5 + a 8 + a 9 = 0. (50) 

These equations do not fix the constants uniquely. We are free to make 
some choices. The term (Foi) 2 in eq.(44) is not present in the original 
Hamiltonian form of the Lagrangian (4), so we choose 

a 3 = 0. (51) 
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We can also arrange for the combination 

1 



^(Ei^ + Ej-Foi (52) 



to appear in as it does in C . This requires (from (50)) 



a2 = c — 2a5 
1 

a 4 = --c + a 5 
4 

ae = gC + a 5 



a 7 = --c - a 5 



a 8 = --c + a 5 
6 



ag = -c 
6 



a = --c - a 5 (53) 



and so 



4 <0 = -4a 1 [-^(F ij ) 2 - i(E,) 2 + E, • F 0t ] (54) 

proportional to the non-ghost part of the original Lagrangian (3). 

Equation (54) does not come from the BRST identities, it just emerges from 
the numerical values of the divergent integrals. It may be a consequence of some 
hidden Lorentz invariance. 

The constants ao,a\,... are still not uniquely fixed. There are two particu- 
larly simple choices. 

(i) Choose ao — with 05 = — f c. Then we find 

11 

ai = _ 12 C 
11 

a 2 = -c 
11 

a 4 = — -c 
6 

a 6 = a 7 = 

5 

as = —c 

a 9 = \c. (55) 
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(ii) The second choice is a\ = with 05 = \c. Then 

11 

6 

a 2 = 

a4 = 
11 
6 
11 

6 
2 

a 8 = --c 

a 9 = -a (56) 

Note that ao has the expected value for coupling constant renormalization. 
The counter-terms in either case are 

A = -^c(F^) 2 - icFtj ■ djAi + ^cgFij • (A, A A,) 

-ic(F 0j ) 2 + ^c(E 4 ) 2 + ^cE 4 • F 0l 

4 4 4 

+ -cgE 4 • (Ai A A ) - -cE, • d A l + -c(u; + d t c*) ■ d t c. (57) 

The counter-terms in a.5, a.6, 07, a§ and ag are involved in a rescaling of the 
fields. Defining 

A^ = (l + a 5 )Aj 
A = (l + a 6 )A 
E m = (1 + «s)E m — a 9 F 0m 

u 'i = (! - a 5)Ui 

u = (1 - a 6 )u 

c' = (1 — 07)0 
K' = (l + a 7 )K 

.9' = (1 + a )g 
c'* = (1 - a 5 )c* 

v' = (l-a 8 )v, (58) 

we have from (48) that 

Co + A = (1 - 4a 1 )£ (</, AJ, A , E', c', c'*, uj, u , K'). (59) 

Note that a6 which determines the renormalization of the Coulomb field Aq 
has the same numerical value as ao- 

We have not calculated the divergences in graphs with four external lines. 
We assume they will be cancelled by the same counter-terms. 
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5 Comments 



We conclude that there is no difficulty to one-loop order in renormalizing the 
Hamiltonian form of the Coulomb gauge. We guess that the renormalization 
would formally go through to higher orders, but then there is the problem men- 
tioned in [1] , [S] , [5] of combining the renormalization of ultra-violet divergences 
with the resolution of energy-divergence ambiguities. 

It is not quite obvious how the renormalization would be formulated if the 
Aq field had been eliminated to give the non-local colour Coulomb potential 
(note the non-zero value of the Aq field renormalization constant ag in (56)). 
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Appendix A 

Here we give as an example the evaluation of the ultra-violet divergent part 
of the graph in Fig. 20. 



xT rz (P)T zv (P)T ru (Q+P)[(-2Q-P) v S uj + (Q-P) u 6 jv + (2P+Q) j S vu }. (Al) 



r (20)abc 
Ojk 



{Q,-Q, 0) = ig' 




PoPk 1 



(p 2 + irj) 2 (q + p) 2 + if] 



Applying the integral 




5 f 

= ^r(-)<7 / dyy{\ - y){(P + yQ) 2 + y(l - y){-q 2 - (A2) 



* Jo 

and power counting to (Al) 




leading to 



V$ )abc (q,-q,0) = ~c 9 r bc qo 5 jk . (A4) 



Appendix B 
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Example of self-energy evaluation H < i^ ab in eq.(ll). Let p, q be internal and 
k external momentum, p — q = k. The sum of two graphs is 

4 T t3 (P)r,(Q) I + {ipo){iqo)]§ab (Bl) 

p z q z 2 

where we have symmetrized the first term in P, Q. the minus sign in the second 
term comes from the opposite order of the f abc factors at the two vertices. Doing 
the po integration by Cauchy, we get 

( 2 ^ 4 ( 2n »i(PT^ (p + Q)[p2 + ° 2 ^ 2PQ]Sab - m 

The last factor (P — Q) 2 is approximately (P ■ K) 2 /P 2 . With this factor, the 
integral is only logarithmically divergent, and to get the divergent part we can 
put Q = P everywhere. We use r iJ (F)T ji (P) = 2. Then we get 

So the divergent part isQ 

UcK 2 6 ab . (54) 
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1 / K iKj 



k 2 + in \ 1 K 2 
j 

'k 2 

E '» Aj ifco /- KmKj 

1 W + i v \ mj k 2 
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K 2 I K m K n 

-2 i :„ \° mn IS2 



k l + a] \ K 



E m A 



iKm 
K 2 



~~* ~lC 2 

Figure 1: Feynman rules for the propagators in the Coulomb gauge. 
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/ c, k d : l \ 



-g 2 f 9ac r b " (V« - W 



c,k 



/>. 



gf^k 



Figure 2: Feynman rules for the vertices in the Coulomb gauge. The arrows 
denote the directions of the momenta. 
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1 




Figure 3: Feynman rules for ghosts and sources in the Coulomb gauge. Doubled 
lines denote ghosts. The black arrows distinguish between ghosts and antighosts. 
Momenta flow into the vertex. 
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Figure 4: The transverse gluon self-energy graphs. 




Figure 5: The A^Aq two-point function. 
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Figure 6: The time-time component of the gluon self-energy. 



Figure 7: The transition between the transverse gluon field and its conjugate 
field Ei. 



S X 



Figure 8: The transition between the Coulomb field Ao and the conjugate field 
E 4 . 



/ N 
J I 



Figure 9: The conjugate field self-energy. 
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Figure 10: Graph contributing to the three-gluon vertex function. 
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Figure 11: There are three graphs in this class with permutations of the vertices. 




Figure 12: Graph representing a class of 6 diagrams. 
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^ ~~ \ y 



Figure 13: There are 3 graphs in this class of diagrams. 




Figure 14: Graph with two external Coulomb lines (there are 3 diagrams in 
this class). 




1G 




A' 



/ 



/ 



/ 



< 



\ 



\ 



NL. 



Figure 17: The graph with two external Coulomb lines and one three-gluon 
vertex. 



Figure 18: Graphs contributing to the(AjA,-Ao) three-point function. 



17 



A 

/ 



Nl 

Figure 19: Graph contributing to the (A^A^Ao) three-point function which 
contains a three-gluon vertex. 
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\ 



Figure 20: The (AiAj A ) graph with a three-gluon vertex. 



v. 



/< 



Figure 21: The (AjAjAo) graph with a four-gluon vertex. 
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S N 

Figure 22: The ghost self-energy. 




Figure 24: The ghost vertex graph with a K source. 
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Figure 25: Graph with external A,;, ghost and anti-ghost lines. 



J! 

Figure 26: Graph with Uq source, E, and c lines. 



Figure 27: Graph with Ui source, Aj and c lines. 



Figure 28: Diagram with v n source, Aq and c lines. 
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I 

\ I 
\ 

Figure 29: Graph contributing to the (A^E^-Ag) vertex function. 



Figure 30: Graph with external gluon, Coulomb and E-field. 
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Figure 31: Graph in the (A,EjAo) vertex function. 
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